See-saw protocol for Doppler cooling of multilevel systems with coherent pulse trains 
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We consider decelerating and cooling an ensemble of multilevel atoms and molecules with a 
coherent train of ultrashort laser pulses. In the frequency domain such trains form frequency combs. 
We propose a novel see-saw scheme of Doppler cooling of multilevel atoms and molecules. In this 
scheme the teeth of the frequency comb are periodically moved in and out of resonance with the 
allowed transitions. The see-saw cooling may be carried out in practice by switching carrier-envelope 
phase offset between predefined values. We study the performance of the see-saw scheme for the 
simplest prototype three-level A system and demonstrate its advantage. For these systems we 
demonstrate a procedure for determining optimal train parameters, and see-saw switching period. 
We also illustrate the performance of the protocol by numerically simulating time evolution of 
velocity distribution. 

PACS numbers: 37.10. De, 37.10.Gh, 42.50.Wk 
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I. INTRODUCTION 

The most widespread method of direct laser cooling of 
atoms is Doppler cooling |l| , in which a sample of atoms 
is irradiated with laser light. Each atom absorbs laser 
photons and then radiates the light in random directions 
returning to the ground state. Repeating this optical 
cycle typically tens of thousands of times can cool the 
atomic sample down to the mK temperatures. Direct 
application of the conventional Doppler cooling to mul- 
tilevel atoms and molecules is challenging: most atoms 
and all molecules have transitions that can radiatively 
branch out to a multitude of other states. Exciting pop- 
ulation from all these lower-energy states requires a large 
number of lasers which makes the conventional scheme 
impractical. 

Direct cooling of only a few special multilevel atomic 
and molecular systems have been demonstrated so far. 
Recently the cooling of multilevel atoms Er and Dy 
were experimentally demonstrated [2|, |3[ ■ Doppler cool- 
ing of a certain class of a molecules with highly diagonal 
Franck-Condon matrixes has been demonstrated The 
prospects for laser cooling of TIF molecules have been re- 
ported 

The enumerated multilevel systems have very unique 
properties; here we address the possibility of cooling gen- 
eral multilevel system. The basic idea is to employ a fre- 
quency comb (FC). Frequency combs are produced by co- 
herent pulse train emitted by mode-locked lasers. Their 
frequency spectrum consists of a regular comb of sharp 
lines, i.e., effectively a large collection of CW lasers. By 
moving the teeth in resonance with individual transitions 
one can selectively transfer populations between the lev- 
els, thereby enabling cooling multilevel systems. 

We start by introducing basic parameters of pulse 
trains and frequency combs. The electric field of the 
train may be expresed as 



E(i) = iE p 22 C0S M + g(t - roT) , (1) 



where e is the polarization vector, E p is the field am- 
plitude, $ m is the phase shift, w c is the carrier laser fre- 
quency and T is the pulse repetition period. $ m = mxtp 
is the accumulated phase with cp being the carrier enve- 
lope phase offset (CEO) between subsequent pulses. The 
pulse envelope g(t) is normalized so that max g(t) = 1. 
The spectral width of a single mode is ~ \/2A~/(NT), 
where N is the number of pulses in the train [6|. For 
an infinite train, N — > oo, the spectral profile of the 
pulse train exhibits a characteristic periodic structure 
composed of a series of regularly spaced laser modes 
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where n is the index of the mode. The frequency interval 
between the neighbouring modes is determined by the 
pulse repetition rate Lo rcp = 2ir/T. The third term, i.e., 
<p/T, shifts the entire FC while keeping intervals between 
teeth constant. 

Our method is based on varying the carrier-envelope 
phase offset ip (see Eq. ([2])). We alter the pulse train in 
such a way that the teeth move in and out of resonance 
with specific transitions. The time-sequence is shown in 
Fig. [T] For illustrative purposes, we consider a four-level 
system with three ground states (|g 1 ),|g 2 ) , |g 3 )) and one 
excited state |ei). Suppose the atom/molecule is initially 
in the \g ± ) state. We start by tuning the comb to the fre- 
quency of the |g x ) — |ei) transition. We Doppler cool the 
atom on this transition. At each radiative emission event 
a part of the population branches out to the |g 2 ) and |g 3 ) 
levels. When the population of the |g 1 ) level is depleted, 
by varying the carrier-envelope phase ip we tune the comb 
to the frequency of the |g 2 ) — |d) transition and cool on 
that transition. Again due to the radiative branching, 
the population will leak to the |g 1 ) and |g 3 ) levels. We 
then retune the comb to the |g 3 ) — |d) transition and 
cool on that transition. And so on. As a result the atom 
or molecule is cooled with a single laser source. We refer 
to this method as the see-saw scheme. 
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FIG. 1: See-saw cooling protocol of multi-level system with 
a tunable frequency comb. 

(a) cp r 

|e) — 10 




FIG. 2: See-saw cooling protocol applied to the A system, (a) 
Shows the energy level frequencies relative to the FC teeth 
positions for ip — (pi. (b) Shows new detunings of the ground 
states when the CEO phase in the subsequent period is set to 
ip = ip2- Alternating the phase between ipi and ipi moves FC 
teeth in and out of resonance with either |gj or |g 2 ). 



II. SEE-SAW DOPPLER COOLING APPLIED 
TO THE A-SYSTEM 

Below we demonstrate how the proposed see-saw cool- 
ing scheme can be applied for cooling the simplest pro- 
totype multilevel system: three-level atoms. We limit 
our consideration to the A configuration, as sketched in 
Fig. [3Jb). The system has two ground states | g x ) and 
|g 2 ), and one excited state |e). The energy gap between 
the two ground states is Ai 2 . The excited state life- 
time r = I/7, where the total radiative decay rate is 
the sum of the decay rates to the two ground states 
7 = 7i+72- The Rabi- frequencies of individual tran- 
sitions are denoted as fii and f^- 



The atomic beam is irradiated by the counter- 
propagating train of laser pulses Eq. ([1} as illustrated in 
Fig.ETa). 

The optical Bloch equations for the density matrix gov- 
erning the evolution of populations p ee , Pg jgj and coher- 
ences p egj , p gi92 (j = 1,2) read 

2 

Pee = -JPee-^2lm[n egj p egj ] , (3) 
3=1 

2 

Peg, = -^Peg, + * ^ ~T~ ^ Pee ^P ~ P 9v93 ) > ( 4 ) 

i=i 

1 

Pg 3 g 3 > = 5 J3'ljPee + ^l g .,P gj e ~ Cl egj Peg^,)- (5) 

The time- and space-dependent Rabi frequency is 

N-l 
rn— 

where Sj — uj c — uj egj , k c — uj c /c and z is the atomic co- 
ordinate. The peak Rabi frequency is fi^ eafc = -j£(e|D • 
e\gj). Eqs. (J31I11E]) were derived using the rotating wave 
approximation. The pulse areas corresponding to indi- 
vidual transitions are: 8i = fl^ eak g(t)dt. The resid- 
ual detunings of transitions \gi) — > |e) from the nearest 
FC teeth, Si can be expressed as 

Si = (fji + 2nn i )/T 1 (7) 

where 

fji = (Si + k c v)T + (p (8) 

is the Doppler-shifted phase offset between the subse- 
quent pulses @, Q , k c is the wave vector and an integer 
m is chosen in order to normalize the detunings to the 
—uj rep /2 < Si < io rep /2 interval. 

Below we assume that individual pulse areas and ra- 
diative decay rates are equal: Q\ = #2> 7i = 72- Di- 
rect population transfer between the two ground states 
is dipole forbidden. Therefore, the second ground state is 
accessible only through the excited state. Additionally, 
the atoms are assumed to be initially prepared in the first 
ground state |g x ). 

The scattering force exerted on atoms by the pulse 
train can be expressed in terms of the coherences as Q 

2 

F z = Hk c ^ lm \Pegj Qe 9] } ■ (9) 
J=l 

In the limiting case when the duration of the pulse is 
much shorter than the repetition period T, one can uti- 
lize the delta-function pulse shape model to describe the 
dynamics of the system Q. In this approximation, the 
fractional momentum kick due to the interaction of an 
atom with the n-th pulse reads 

AjWpr = (Pee)r ~ , (10) 
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(a) Pulse train 




FIG. 3: (a) An atomic beam of A-type systems interacting 
with a counter propagating pulse train, (b) Energy level di- 
agram of A system, (c) Frequency domain picture of comb 
interacting with A system. 

where p r — hk c is the recoil momentum, and (p"e); r are 
the populations of the exited state right before and after 
the n-th pulse. 

The average scattering force may be determined as 

F- = (ID 

where Ap is the total momentum imparted by the pulse 
train, T is its repetition period and N is the number of 
pulses in the train. Maximum value of F avg = p r /T is 
attained when each pulse transfers full recoil momentum 
to the atom. 

We suppose that initially the central FC tooth is reso- 
nant with the |e) — > Ig-J transition. Then the frequency 
of the |e) — > |g 2 ) transition will be detuned with respect 
to the central tooth by — A12, where A12 is the energy 
gap between the two ground states. The nearest to the 
uj e g 2 tooth will be detuned from it by ?7io<i(Ai2, lu t 

ep) if 

mod(Ai2, u>rep) < u rep /2 and mod(Ai2, u> rep ) — to r ep if 
mod( A12, uJrep) > u rep /2. In the first case it is blue de- 
tuned and in the second case it is red detuned. The pulse 
repetition frequency is chosen so that mod( A12, io rep ) ^ 
0, that is the two photon resonance condition is avoided. 
In case when mod{ A12, to rep ) = the system evolves into 
the "dark" superposition of the two ground states which 



is transparent to the pulse train [l(| and therefore does 
not exert the scattering force. The population of the sys- 
tem is initially in the first ground state |g 1 ). If the atoms 
are irradiated by the pulse train for some time, the popu- 
lation will migrate from the first, strongly driven, ground 
state \gi) to the second, less coupled, ground state |g 2 ). 
Then according to the see-saw protocol, the phase offset 
can be switched tp = tp± = k c v mp T to <p = 952, where 
9i = k c v mp T - mod(Ai2,uj rep )T if mod( A 12 , ui rep ) < 
uj rep /2 and cp 2 — k c v mp T + uj rep — mod( A12, u> rep )T if 
mod( A12, uJrep) > Wrep/2 ( where v mp is usually the cen- 
ter of velocity distribution), so that the second transition 
becomes resonant with one of the FC modes. 

To evaluate the time evolution of the density matrix 
we numerically integrated the optical Bloch equations (J3J 
21 [5]). Fig. f4][a) shows the time evolution of population 
distribution over the first 100 lifetimes. One can see that 
as the system approaches the quasi-steady state regime, 
the average population becomes redistributed mostly be- 
tween the two ground states and the post pulse excited 
state population saturates without further increase. To 
interrupt the saturation a change in the parameters is 
needed. This is accomplished by changing the phase 
from tp = ipi to tp — if2- Such switching is illustrated 
in Fig. [lb). 

Fig. |4jb) shows the mechanical momentum accumu- 
lated by the atoms. It is clear that at the end of each 
50-pulse period, the momentum accumulation saturates 
to a constant value. This means that the rate of de- 
celeration slows down. Once the phase is set to tp — <p% 
after the 50 tft- pulse, a sudden increase in the accumulated 
momentum occurs which improves the cooling efficiency. 
Continuous switching of the phase between ip — ipi and 
ip = (pa every 50 th pulse sustains more efficient momen- 
tum transfer from the pulse train to the atomic beam. 
Fig. [5] compares the excited state dynamics for the non 
see-saw and the see-saw schemes, switching the phase be- 
tween ipi to <f2 or vice versa every 50*' 1 pulse (i.e., red 
curve) . The average excited state population is increased 
compared to the case when the phase is fixed (i.e., black 
curve) . 

Hereafter in this paper, we denote N sw as the length of 
sub-train for which we keep the phase constant at either 
values ip\ and ip 2 - For example N sw — 50 in Fig. |4j 



III. COOLING A ID BEAM OF A SYSTEMS 

In this section we apply the proposed see-saw cooling 
method for cooling a thermal ID beam of three-level A- 
type atoms. First we find the optimal parameters maxi- 
mizing the scattering force exerted on the atoms. Second 
we study the evolution of the ensemble velocity distribu- 
tion during cooling with these parameters. 
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FIG. 4: (Color online) Time-dependence of A-system popu- 
lations (panel (a)) and fractional momentum kick (panel (b)) 
in the see-saw protocol with N 3W = 50. These figures were 
computed by integrating OBEs, Eq. @ with the following 
parameters: ground states splitting A = 2n x 300 GHz, spon- 
taneous decays of 71 = 72 = 2n x 10 MHz, repetition period 
is T = 0.5/7 e ns i an d pulse areas of 9i = #2 = 7r/10. 



A. Maximum scattering force 

The atoms under consideration have the A-like config- 
uration of Fig. |3l Specific atomic values are taken to be 
as follows: ground state splitting is A = 2ir x 300 GHz. 
Spontaneous decay rate to the ground sates are equal, 
i.e., 71 = 72 = 2tt x 10 MHz, leading to the excited state 
lifetime of 8 ns. Pulse areas are taken to be equal as 
well, 6*i = 62- We optimize the FC parameters for the 
two different values of the pulse area 6*1,2 = tt/10 and 
tt/2. 

All parameters may be categorized into three groups 

• Atomic parameters: A 12, 71 and 72. 

• Fixed laser field parameters: T and laser intensity. 

• Variable laser-field parameter: ipi, tf2 and N aw . 

The description of the optimization process is summa- 
rized in the following text. Parameters A12, 71 and 72 
remain fixed for the specific A system. A mode-locked 
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FIG. 5: (Color online) Excited state population as a function 
of time. The black curve is for a fixed tp. The red curve is for 
see-saw scheme with N 3W = 50. The jump in the population 
is apparent once the phase is switched from ip = to <p = 
— mod(Ai2, 0Jre P )T. These plots are obtained for the following 
parameters: ground states splitting A12 = 2n x 300 GHz, 
spontaneous decays of 71 = 72 = 2n x 10 MHz, repetition 
period is T = 1.00017 ns, and pulse areas 9\ — 62 = tt/10. 



laser source is assumed to have fixed repetition period 
T and intensity. The intensity of the laser beam along 
with the transition dipole moments determine the Rabi 
frequencies f2i,2 and pulse areas 6*1,2- An optical element 
may be placed at the output of the laser source to control 
the CEO phases ipi and ip2, and the see-saw switching 
period N sw . 

We carry out an exhaustive search for the optimal val- 
ues of T and N sw corresponding to the maximum possible 
average scattering force. The maximum and minimum 
pulse repetition period values are taken as T m i n = 1 
ns and T max = 10 ns. Maximum and minimum num- 
ber of pulses between the switching are N™ n = 1 and 
N™ x = 60. As an example, we obtained the optimal rep- 
etition periods and the optimal see-saw switching periods 
for the two pulse areas 6*1,2 = tt/10 and 6*1,2 = tt/2. For 
6*1,2 = tt/10, optimal pulse repetition period was found 
to be T opt = 1.0035 ns and NX = 1 - For #1.2 = tt/2, 
T *" = 1.0083 ns and NX = !• 



B. Evolution of the velocity distribution 

We consider a ID atomic beam with the initial velocity 
distribution 



9 v 3 3u 2 
f(v,t = 0) = -_exp(-— 5-), 



(12) 



where most probable value of the velocity v mp is set to 
300 m/s. 

The carrier laser frequency is resonant with the atomic 
transition frequency, lo c = w egi . The values of the phase 
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FIG. 6: Right axis is the average scattering force F avg plotted 
for the Doppler shifted velocities v in the atomic beam. The 
left axis is the average momentum kick per pulse Ap/(Np r ) 
imparted to the atoms in the same atomic beam. 



offsets ipi, ip 2 are chosen with account of the Doppler 
shift of the laser frequency in the reference frame moving 
with the center of velocity distribution, ipi — k c v mp T, 
ip 2 = k c v mp T - mod(A 12 , io rep )T (if mod(Ai 2 , uj rep ) < 
uj rep /2), (p 2 = k c v mp T + u rep T - mod( A 12 ,u> rep )T (if 
mod( A12, uirep) > oj rep /2). Initially ipi = k c v mp T and as 
the atoms slow down <p\ follows the v mp (t). 



Fig. |5] shows the fractional momentum kick per pulse 
Ap/(Np r ) versus atomic velocity for the pulse area 



tt/10, T 



1.0035 ns and N s . 
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optimal parameters correspond to the teeth separation 
of uj rep ~ 2n x 1 GHz in the frequency space and 
v = A c /T ~ 587 m/s in the velocity space. 
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The FC teeth are initially resonant with the atomic 
transition frequency w e9l for the atoms with velocities 
Vjjn = v mp + n\ c /T. Continuously as the atoms slow 
down the phases tpi, ip 2 have to be shifted by the amount 
—k c Av mp T (where Av mp is the decrease in the value of 
the central velocity due to the slowing of the beam) to 
compensate for the change in phases fjj, Eq.©. This 
leads to the formation of pronounced peaks in the ve- 
locity dependence of the scattering force, separated by 
A c /T ~ 587 m/s. Fig. [6] shows one of such peaks for 
n = or vdq — v mp — 300 m/s. 



velocity (m/s) 



FIG. 7: The plot shows the two velocity dependent peaks of 
fractional momentum kick that were not resolved in Fig. [6] 
due to overlapping of frequency comb teeth. The plot is for 
pulse areas 0i t 2 = 7r/100 otherwise the parameters are the 
same as Fig. [6] 



The velocity dependence of the scattering force can be 
explained on the basis of the Fourier transform of the 
see-saw pulse train EqflJ The Fourier transform (the 
derivation is given in the appendix) of the electric field 
(P) reads 



J 



= F/H sin( ^ (l) xV^- 2Br 



shi(±Ci) 



sin^V^) 
sm(K2) 



2nn , 



(13) 



r 



where Q = cj>i + (u)—u) c )T, ( = (0i + 2 ) + 2(w-w c )T, and 
F{ 2 (iu) are the Fourier transforms of the pulse envelopes. 



the frequency of corresponding atomic transition for a 
given atomic velocity. 



The arguments of the two (5-functions determine po- 
sitions of the FC teeth: ( — j^- — or uj n — uj c 4 

— . In the atomic reference frame these val 
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ues are Doppler shifted: u^-a^Tu-cT ytn 2T~ 

When u° = w eo ., then the n-th tooth is resonant with 



Now we return to the interpretation of the velocity 
peaks in Fig. [7l The carrier laser frequency is res- 
onant with the frequency of the first atomic transi- 
tion uj c — uj ea , . Then for a set of velocities nJJ = 

/k c the FC teeth will be resonant with 
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FIG. 8: Evolution of the velocity distribution of a ID beam 
of A-type atoms interacting with a train of ultrashort laser 
pulses. Curve (a) displays the initial velocity distribution of 
the atomic beam, curve (b) shows velocity distribution after 
interaction with 25,000 pulses, curve (c) shows velocity distri- 
bution evolution after deceleration by 100,000 pulses. These 
were plotted for 9i, 2 = 7r/10, T = 1.0035 ns and N sw = 1. 
The remaining parameters are the same as in Fig. [5] 



the frequency of atomic transition u egi . 
time for the velocities v b „ = (iiffo _ » / 
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another transition. The difference between the velocities 
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< 0.02 




100 200 300 400 500 600 

velocity (m/s) 



FIG. 9: Fractional momentum kick as a function of the 
Doppler shifted atomic velocity for repetition period T = 
1.00083 ns, switching period of N 3W = 1, and pulse areas 
#1,2 = it/2. Otherwise the parameters are the same as in 
Fig. El 



pulse transferred to the atomic ensemble for single pulse 
areas #1,2 = tt/2. The optimal parameters for which the 
scattering force is maximized are: T — 1.00083 ns and 
N s w = 1. 

For the large pulse areas 9\ y 2 ~ tt/v^ the teeth in 
the scattering force spectral profile are essentially power- 
broaden. Low contrast of the maxima in velocity depen- 
dence of the scattering force leads to the slower compres- 
sion of the velocity distribution. As a result the ensemble 
can be decelerated, but not cooled. 



< - < = - nx) + A 12 \ /k c . (14) 

For N sw — 1 the minimum distance between the two 
resonant velocities equals to mod(Ai2, u rep )/k c . If the 
pulse area is small (6i -C 1) and the lifetime of the excited 
state is small compared to the pulse repetition period, 
then the two peaks corresponding to the nearest velocities 
and v„ 2 can be resolved, Fig. [7] 

Figure [3] shows time evolution of the atomic velocity 
distribution for various number of pulses. This is a nu- 
merical simulation for a ID atomic beam with the initial 
velocity distribution given by Eq. (112p . We present two 
computed snapshots of velocity distribution for 25,000 
and 100,000 pulses. For the interaction with 25,000 
pulses, there is a deceleration for a major part of the 
distribution accompanied with velocity compression. As 
the number of pulses interacting with atoms increases, 
the continuous deceleration and the decrease of the width 
of velocity distribution is becomes more pronounced, as 
shown in Fig. [8] curve (c). 

As the pulse area is increased, the scattering force 
ceases to become velocity selective. This is illustrated 
in Fig. IHl where we present the average momentum per 



IV. CONCLUSION 

We have proposed the see-saw protocol for decelerat- 
ing and cooling multilevel systems by trains of ultrashort 
laser pulses. We have demonstrated the efficiency of the 
proposed method in cooling the ensemble of three-level 
A-type atoms. The method is based on periodic interrup- 
tion of the saturation regime in the system by switching 
the phase offset between subsequent pulses. As a result 
the sustained population transfer to the excited state oc- 
curs and nonzero scattering force is exerted on the atoms. 
The see-saw scattering force is velocity-dependent and 
mimics the Fourier transform image of the see-saw pulse 
train. At large single pulse areas = \J d\ + Q\ ~ 7r the 
scattering force becomes largely velocity insensitive. The 
same effect happens when the pulse repetition period is 
much longer then the excited state lifetime. The FC pa- 
rameters, such as intensity, pulse repetition period and 
the number of pulses between the phase switchings can 
be optimized in order to maximize the scattering force. 
Additionally the phase offsets has to be adjusted during 
the deceleration process, while the position of the center 
of velocity distribution moves toward the lower velocities. 

We illustrated the see-saw cooling by computing time 
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evolution of velocity distribution for a thermal beam of 
three-level A-type atoms. The center of velocity distri- 
bution moves toward lower velocities, i.e., the beam is 
slowed down. Moreover the width of the distribution is 
decreased, velocity distribution is compressed, i.e., the 
ensemble in addition to being slowed is cooled. 
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Appendix: Fourier Transform of a Pulse Train with 
see-saw CEO phase switching 

Consider a train of identical pulses emitted from an 
ideal mode-locked laser at equal time intervals. The rep- 



£{t) 
£ x {t) 



where g(t) is the envelope function. The two terms in 
the preceding equation differ by the time shift of N SW T 
and their accumulative phases. We compute the Fourier 



Here F*(cj) is the Fourier transform of the first pulse of 
the first term in Eq. (|A.l|) and F^(u>) is is the Fourier 
transform of the first pulse of the first term in Eq. (IA.2I) . 
For example, in case of a Gaussian envelope, g(t) — 



etition period between the pulses is fixed. We consider 
the case when the pulses accumulates a CEO phase of ifi 
and ip2 every N sw pulses alternatively. Fig. 1101 illustrates 
such pulse train in the time domain. The snapshot of the 
pulse train in Fig. 1101 shows that the train may be subdi- 
vided into N shorter trains. Each shorter train has 2N SW 
pulses. The first N sw pulses accumulate a CEO phase of 
<pi every pulse. The subsequent N sw pulses accumulate 
a CEO phase of <f2 every pulse. Therefore, the phase of 
a given pulse equals to the total accumulated phases for 
all the previous individual pulses. The electric field at a 
fixed position for a pulse train of Fig. [TU] is 



(A.l) 
(A.2) 

I 

transform of Eq. ([A~T|) and Eg. dA~2]) as £(uj) = T[£] = 
—jk= f °!L £(t)e lu ' t dt. The result with n = uj — uj c reads 

V ^TT J 00 



(A.3) 



In Ea. (|A.3|l we ignored the Fourier transform of the c.c. 
terms in Eq. (|A.l|) and Eq. (|A.2[) as we are interested only 
in the frequency comb that is going to be built around uj c . 
By treating each sum in Eq. (|A.3[) as a geometric series, 
we obtain 



=£ 1 {t)+£ 2 {t), 

1 N-IN SW -1 

m— n— 
N-lN, m -l 

.±£q (g(t)e~ i " ct+i ( n+mN "" +1 ^ V2+l ( mN "" +Nsu '~ lS>iP1 +C.C.). 



+ c.c), 



m—0 n— 



JV-1 N sw -l 

£(u) =[Ff(u) ^2 ^2 e l ^ n+2mNs ^ vT+{n+mNsm)ipi+mNs, " V2) 

m—0 n— 
N-l N sw -l 

m—0 n— 

i 
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FIG. 10: The static picture of a pulse train that accumulates phases y>i and ipi alternatively each N 3m pulses. 



sm(±N sw (iT) sm(±NN sw {T) sin (±A^C 2 !T) sin {\N N sw (T) 

+ t 2 (u)- 



sin^Ci? 1 ) sm{±N sw (T) 



sin(iC2T) sin(±N sw CT) 



(A.4) 



where Q = fa/T + (lu - lo c ) and £ = (<£i + <fe) /T + 2(w - 
cj c ). In the limit when TV — > 00, Eq. (|A.4p simplifies to 



g H=F 1+H sin( ^; c f ) 



1 v atu J 



shijlAUCgT) 
sin(iC2T) 



2n7r 



(A.5) 



[1] T. Hansch and A. Schawlow, Opt. Comm. 13, 68 (1975). 
[2] S. H. Youn, M. Lu, and B. L. Lev, Phys. Rev. A 82, 
043403 (2010). 

[3] M. Lu, N. Q. Burdick, and B. L. Lev, Phys. Rev. Lett. 

108, 215301 (2012). 
[4] E. S. Shuman, J. F. Barry, D. R. Glenn, and D. DeMille, 

Phys. Rev. Lett. 103, 223001 (2009). 
[5] L. R. Hunter, S. K. Peck, A. S. Greenspon, S. S. Alam, 

and D. DeMille, Phys. Rev. A 85, 012511 (2012). 
[6] T. Haensch, Metrology and fundamental constants (IOS 



Press, 2007). 

[7] E. Ilinova, M. Ahmad, and A. Derevianko, Phys. Rev. A 

84, 033421 (2011). 
[8] E. Ilinova and A. Derevianko, 

http://arxiv.org/abs/1203.1963 (2012). 
[9] E. Ilinova, M. Ahmad, and A. Derevianko, Phys. Rev. A 

84, 033421 (2011). 
[10] L. E. E. de Araujo, Phys. Rev. A 69, 013408 (2004). 



